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Abstract 

The thermodynamics of a lattice regularized asymmetric Anderson impurity in a correlated host is 
obtained by an exact solution. The crossover from the Anderson- to the Kondo-regime is studied, thus 
making contact with predictions by scaling theory. On the basis of the exact solution, the transition 
to universal Kondo behavior is shown to be realized by a graduate separation of the energy scales of 
spin and charge excitations. 

1 Introduction 

Homogeneous correlated one-dimensional electron systems have been studied intensively in recent years 
and exact results have been obtained by as different techniques as bosonization [1,2,14], conformal 
field theory (CFT) [11], and the algebraic Bethe Ansatz (ABA) [24]. Whereas bosonization and CFT 
deal with the ground state and low-temperature properties of the continuum models, the Quantum 
Transfer Matrix (QTM) technique has been developed [22] for the calculation of the thermodynamics 
of lattice models. The QTM combines the Bethe Ansatz solution with the quantum inverse scattering 
method and a Trotter-Suzuki mapping of the partition function on a classical two-dimensional lattice. 
In this approach, the free energy of the model is encoded in a set of finitely many non-linear integral 
equations (NLIE). All these methods have been applied for example to the general spin-1/2 Heisenberg 
chain [22], the t J-model [19] and the Hubbard model [21,23]. At low temperatures, i.e. in the conformal 
limit, the different methods yield identical results. 

The equivalence of continuum and lattice methods is much less investigated in the case of inho- 
mogeneous systems, where a localized impurity interacts with the host by s-wave scattering, so that 
the model is reduced to one dimension. First, the Kondo and the Anderson models have been solved 
exactly by the coordinate Bethe Ansatz [30], in this case starting directly from the continuum model, 
i.e., a point-like impurity coupled to a free host. This approach has been extended to calculate ther- 
modynamic properties by applying the Thermodynamic Bethe Ansatz (TBA), which yields the free 
energy in terms of infinitely many coupled NLIE. In the case of the spin-1/2 single channel Kondo 
model, the presence of a stable infrared fixed point is the starting point for CFT-calculations [1,2]. 
These confirm results of the Bethe- Ansatz solution and go beyond. 

Later, several schemes of including an impurity in the correlated host on a lattice accessible by 
the QISM have been found [5,6,8]. However, the connection between the continuum and the lattice 
models in the conformal limit with a generic Kondo-like interaction has not been established until 
very recently [9]. In [9], a lattice regularized asymmetric Anderson impurity in a correlated host has 
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been obtained from QISM, and the free energy has been calculated exactly by the QTM-method. As 
a special case, the Kondo limit has been performed explicitly, yielding a conformally invariant bulk of 
free fermions interacting with a magnetic impurity through spin exchange. The results in the Kondo 
limit are identical with and extend those from the BA-solution, but they are derived from a finite 
number of NLIE. 

Here, we study the properties of the model proposed in [9] in the interesting regime of crossover 
from Anderson- to Kondo-like properties. The model consists of a three-state impurity permitting 
single and zero occupation in a correlated fermionic host. The impurity interacts with its nearest host 
neighbor lattice sites by charge fluctuation and spin exchange. By employing the Poor Man's scaling 
approach proposed by Anderson [4] in the framework of the Kondo model and applied subsequently 
to the generic asymmetric Anderson model [17,18], three different regimes at low temperatures are 
qualitatively distinguished, namely the mixed valence regime (MV), the strong coupling regime (SC) 
and the local moment regime (LM). This qualitative description is confirmed quantitatively on the 
basis of the exact solution. Especially, the universal strong coupling limit is attained by energetically 
decoupling charge and spin excitations. 

The key ingredient of our approach is the one-parametric family of four-dimensional irreps of 
gl(2|l), from which a gl(2|l)-symmetric transfer matrix is constructed, yielding the impurity model 
in the Hamiltonian limit. We use different representations labeled by different interaction parameters 
in the host and impurity spaces. The motivations to consider gl(2|l) are twofold: On one hand, the 
limit of free fermions is contained as a special case and is given by algebraic functions (rather than 
by trigonometric functions like in the XX Z model). On the other hand, su(2) is one of the even 
subalgebras of gl(2|l), and one expects that gl(2|l) symmetry can be reduced to su(2) symmetry in a 
certain limit of the parameters. This is indeed true for special parameter choices and constitutes the 
Kondo limit. 

This article is organized as follows. In the next section the model is defined, including a qualitative 
analysis in the Poor Man's scaling approach. The third section contains the quantitative study of the 
different regimes on the basis of the exact solution. The results are summarized in a phase diagram. 
The fourth section focuses on the Luttinger-liquid behavior of the host. In the last section, conclusions 
are drawn. 

2 Definition of the model and qualitative analysis 

The Hamiltonian is defined as the logarithmic derivative of an inhomogeneous transfer matrix given 
in appendix E It is the sum of the host and impurity contributions [9], 

H = Hh + PH,F + He^ (1) 
Hh = -J2i'^D{cosk-l)+fi)cl^Cka + Hint (2) 

k,cr 

Hi = Eq- {Eq + fl)nd + ^ [Vfc4^dCT + /i-C.j + ^ Jk,p (^4_^Cfc5^dld^ + cl aCkanda^ (3) 

k,a- 



k,p,a 
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where a = —a. Interactions in the host are allowed which are parameterized by a representation 
theoretic parameter a such that Hint = C'(l/a) [7]; the limit a ^ oo corresponds to free fermions. 
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The limit a — > yields the tJ-model, which can be seen by a canonical transformation [7]. In 
the impurity contribution is projected onto non-doubly occupied states by the projection operator P. 
The energy of the non-occupied impurity state in Q is i?o = — 4L'a^/(4uQ + a^). In the limit a ^ 1, 
the coupling constants are given by 

Jk,p = TTr^^^fi^ofe'^'-e-^pUfe''= + e"^f)V (5) 



2{ul + ay 4) 

Additionally to [9], where we focussed on the most important terms even in uq, the leading terms 
uneven in uq are included here. 

The model describes an asymmetric Anderson model in a correlated host, where the doubly- 
occupied impurity state is dynamically decoupled and energetically suppressed. Interactions are tuned 
by a,uo, the bandwidth D and the chemical potential fi. The homogeneous system, which consists of 
Hh alone, was studied in [20,27]. To obtain a qualitative understanding of the model, note that for the 
impurity model with uq = 0, the thermodynamic quantities depend only on T/D for fixed n/D. The 
meaning of different choices of for the finite-temperature behavior is revealed by studying the 
scaling properties following the Poor Man's approach, [4,17,18]. It consists of investigating the leading 
renormalization of the coupling constants by the reduction oi D to D — \6D\ =: D^^^^\ integrating 
out states with energies |e| > D^^^^\ This renormalization procedure is done in the limit T ^ D, at 
sufficiently low temperatures compared to the bandwidth. The resonance width A{lo) is defined by 

^•^^ ^ ~ ^^^^'^ ' 

fc = — TT 

where = 2D{cosk — 1) + /i is the bare one-particle dispersion relation of the host, eq. ((SJ, and 
rik is the occupation of the energy level e^. The resonance width has a maximum at the Fermi level, 
Aq = A, and decays to zero at the band edges. Let us consider the scaling properties of A, and the 
energy difference between single and non-occupation, 

£d ■= - Eq. (6) 

By performing a canonical transformation, the hybridization terms in ^ are eliminated to leading 
order, see [9] for details. Let us first concentrate on the relations between 6D and 6ed, (5A (constant 
prefactors are omitted): 

5ed ~ -A6lnD 

A • Aon 
SA ~ — -^^SlnD 

From the last equation, A remains unrenormalized in the limit D — > oo. The first equation can then 
be integrated to give e^^^'^^ + Aln D^^'^'f') = e^, where is a scaling invariant and e^^'^^^ := ed{D^^ff^). 
This means that e"^^^^ grows logarithmically if D^'^^f^ is reduced. Depending on the relation of e^^^'^ 
to D^^ff \ one distinguishes between three cases: 

i) ef'^^^ e [-D^^ff\D^''ff'^]: Mixed valence (MV) 
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(eff) 



> L>("//): Strong coupling (SC) 



iii) Completely filled band such that D^^ff'> = or D^^^^'^ is not defined: Local moment (LM). 

The first two cases are realized by ^ sufficiently large to ensure a finite filling of the band but fi < jic, 
where the critical potential is calculated below from the exact solution as Hc = 4Z)(a + l)/a. In the 
MV case, both /i and D are finite, so that s^^"^^^ is situated in the band or near the band edges with 
the impurity orbital retaining its triplet degeneracy. The signature in the susceptibility is a crossover 
from Tx{T) = 1/6 at temperatures T > A to the Fermi liquid property xi^ = 0) =const., with 
a non-universal constant value. It originates in a non-integral impurity occupation number at low 
temperatures; that is why this case is referred to as the "mixed valence" (MV) regime. 

The conditions of the second case are met for /i < /ic in the limit /x ^ oo at fixed /x/D. Single 
occupation of the impurity is enforced and ^^^^^^ moves outside the band upon scaling, so that only 
virtual charge fluctuations between the impurity and the host occur. This crossover takes place at 
an effective bandwidth D* given by D* := —e'^^^'' = D^^ff\ i.e. D* — Aln_D* = —e*^. The exact 
solution presented below shows that D* decreases with decreasing D/ ^. Since here, we are dealing with 
temperatures T D, where D is the initial bandwidth, one can neglect the momentum dependence 
of Jk^p and take the momenta at the Fermi points Jk,p = Jkp,pp- The linear order in uq contributes 
only if pp = kp, so that we define Jkp,kF —'■ J- Now the Poor Man's approach is applied to the spin 
exchange terms ©, which again can be done by a canonical transformation eliminating these terms 
to the leading order. This results in 

6J -J^dlnD, (7) 

from which a new scaling invariant 

Tk := De-y' = D^^ff^e-'^'''"'' (8) 

is deduced, where J is given by @. In the limit \uq\ ^ 1, the sign of the spin exchange depends 
on both the sign of kp and uq. For no < 0, spin-exchange is antiferromagnetic at the right Fermi 
point kp > and vice versa. From Q it follows that antiferromagnetic exchange dominates over 
the ferromagnetic term in the scaling procedure, so that from ©, the antiferromagnetic exchange 
constant is given by J ~ | sin A;i?|/|iio| for \uq\ ^ 1. Comparison with Q leads to \uq\ ~ InD. The 
susceptibility is expected to show a crossover from Tx{T) = 1/6 to 1/6 < Tx{T) < 1/4, ending at 
x(0) = 1/{2ttTk), which is the universal value in the strong coupling (SC) regime [31]. 

A further reduction of the initial bandwidth D, or an increase of the chemical potential at fixed D 
such that fj, > fic moves D* and drives the system into the third case. Then at low temperatures, 
the chemical potential lies outside the band. No interactions between the impurity and the host are 
possible, so Tx{T) changes from 1/6 to 1/4 from high to low temperatures. This is the local moment 
(LM) regime. 

Summarizing the results from scaling arguments, one finds that uq plays a crucial role in the 
crossover to the strong-coupling limit and is determined by the renormalization of the antiferromag- 
netic spin-coupling constant J. The additional no-dependence of A, see @, is neglected in the Poor 
Man's scaling; however, it has to be included in order to render the model integrable. 

Let us briefly compare our model with the generic asymmetric Anderson model [3], 

Hi = ed'^ Ud,^ + U ndiUdi + ^ Vkcl^d^ + h.c. (9) 
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where U ^ |£d|,A such that double occupation of the impurity site is excluded. Contrary to @, 
the model does not include backscattering [26]; it contains a bare spin-exchange term, which is 
generated in Q by a canonical transformation [29] in the limit U/A ^ 1. These modifications of 
Q are necessary in order to make our model integrable on a lattice, but do not alter the physics: 
The model © also exhibits the three regimes discussed above in the appropriate asymmetric limit 
U^ed,A [17], with the further possibility to tune in @ such that the impurity is completely 
depleted at low temperatures. With respect to the three regimes discussed in this work, and Q 
are equivalent. 



3 Exact solution 

The rather qualitative picture of the preceding section is made more precise by the exact solution 
yielding the free energy of the host per lattice site fh and the impurity fi in terms of suitably chosen 
auxiliary functions 53 = l + b,*B = l + b,^=l + c, cf. appendix EI 



$ * In ■ 



$ * In — 
c 



lnb{v) - 
lnb{v) - 

lnc{v) = -pD{a + l)^a{v) + PfJ, 

+ [^> * ln^]{v - i/2) - * ln^]{v + i/2), 
where the integration kernels and driving term 



{v + i/2)+f3h/2 + [k*ln^]{v) - [k*hi^]{v + i) 
{v - i/2) - ph/2 + [k* ln^]{v) -[k* ln^](t; - i) 

e:] 

A; * In — {v) 
c 



1 



2 sinh TTV 

f oo 



1 



k{v) = — 



p-|fct/2 



a 



2tt 2coshA;/2 



f 2 + q2/4 



e^'^'^dk 



(10) 
(11) 

(12) 

(13) 



have been defined. The convolutions are defined as [/ * g]{v) := f^^f{v — w)g{w)dw. With the 
solution of (fTn)) - H12() one obtains 



fi 

fh 



-2fi + 2D-T 



In - * k + Dk 
c 



(^^o) 



2tt 



[In 5353] (t;) 
cosh(no — v) 



dv 



(14) 



^fQ, * In - 
c 



(0) dv. 



In this section we discuss the three regimes distinguished above. The NLIE are therefore solved 
analytically near complete filling at zero temperature and numerically for arbitrary parameters. 

At low temperatures the driving term of Inc, eq. ()12() . gets arbitrarily large and determines the 
location of zeroes of In c. It is convenient to define the quantities fl, /ic by 



4D{a + l)/a. 



(15) 



It is shown below that: 



5 



i) No zeroes of In c occur for fi > ^Uc- This case corresponds to complete filling of the band, where 
the impurity is in the LM regime. 

ii) In c acquires two zeroes it In for fj, < fic- For finite 

\uo\ = (lnZ))/7r, (16) 

I no I ^ Ac, the MV regime is realized. From (|16|) . justified rigorously in section [3.31 and the 
D-dependence of A^ calculated below ((T7j) . one identifies ~ \uo\ and D^^^^^^ ~ A^. 

iii) If in the MV regime uq in is chosen such that |no| ~ A^ (corresponding to d'^^^^^ = D* 
defined in the previous section), the crossover to the SC regime takes place. This regime is 
completely realized in the limit — > oo at fixed n/D. 



3.1 LM regime 

For /2 < at low temperatures, InC^ = Inc. Then the eqs. (jlUj) . are trivially solved by In b = 
Ph = — Inb, such that the free energy p4)) is essentially that of an uncoupled local moment. Thus 
XiiT ^ 0) = 1/(4T) diverges at lowest temperatures. At high temperatures, Inb, Inb, Inc approach 
their asymptotic values [9], such that the susceptibility reaches the value xi^ — > oo) = 1/(6T). The 
crossover between these two extremes is shown in fig. with the values D = 1, 30 at /U = 300. 



3.2 MV regime 

To gain analytical insight at T = 0, scaled auxiliary functions and their zeroes are defined: 

et, = rinb ; e, = Tlnc 
eb(±Ab) = ; e,{±A,) = 

Since in the low-temperature limit Inb — > — oo at finite h, the function IniB vanishes. Following [10], 
one can apply Wiener-Hopf-techniques to the linearized NLIE (appendix^ at low temperatures for 
1/2, /i| <^ 1. Performing analogous calculations as in [10], one can show that: 



3/2 

e,{v) = -D{a + l)^a{v) + f^-—cosli{Trv) -j^k{v)jl^''^ 

^ Q-K{D{a + l)Y'^ 

+ ._,,T''\..r,, n + 0{a-,^l\-,^l-) (18) 



8^2 (D(a + 1)) 
^ := / e™ec(w^)du; 
64a + 1 



3/2 

/i3/2 + 0f^5/2^ (19) 



3 a3 



Kl + O (A^) . (20) 



Here a := ^p- and we defined a quantity 7 whose meaning will become clear later. The above 
expressions are valid for small /2(> 0) and /i, i.e. near complete filling and small magnetization. Note 
that A^ is not defined for /i < 0, and neither is 7. 
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Inserting the above results into the expression for the free energy 1)14(1 . (|15() one obtains: 



„ 7cosh(7rno) a^''^k(uo) _.:>/2 a^''^jk(uo) 

ti = —fJ- a TTTM H TTTTflM 

^ 27r(D(a + l))^/2 Svr^ (^^la + 1)) ^ 

n.(/. = 0) = l--4^^^^^^/^ = l-^n. (21) 
cosh(7rno) ^n. 



xih = 0) 



7r7 67r7 

2D{a + l) _ 2_ / g \ 3/2 , ^ 



1/2 / \ 1/2 



a 67r \D{a + 1)/ vr^ \iJ(Q! + 1)/ 

nh{h = 0) = 2-- " =: 2 - 5n;, (22) 

vr yL/^a + 1)/ 

X.(/i = 0) = ^. (23) 

At constant fi, the impurity occupation and susceptibihty are tuned by the parameter uq: For finite 
uq, the impurity site is not fuhy occupied, and neither is the host since /x > 0. The susceptibihties are 
expressed by the deviation from complete occupation, which justifies the term "mixed valence". For 
uq ^ 00, single occupation on the impurity is enforced, whereas the host is away from complete filling. 
Especially, for uq = (lnD)/7r, the impurity susceptibility is independent of D. It is shown below that 
this choice of uq leads to the Kondo regime for all temperatures. 

The high temperature asymptotics are Xii^ 00) = 1/{QT) and Xh{T — > 00) = 1/(8T), corre- 
sponding to three and four available states, respectively. 

The crossover of Xii^) from low to high temperatures at arbitrary values of p, is shown in figs. ^ 
and El 

3.3 SC regime 

By increasing D at fixed n/D and with uq = (lnD)/7r, the crossover from the MV to the SC regime 
is performed. In terms of the zeroes itAc of Inc, this is the case when A,; <C (lnL))/7r. Then by 
substituting v = —x + (lnD)/7r, the convolution in the driving terms in (fTU)) . is approximated by 

[In ^*^]{v) = -Pe-^'^TK , (24) 

where 

/oo 
e™ln[e:/c](u;)du; =: K. (25) 
-00 

The formula for Tx in [9] is an approximation of ()25() . From ()25|). it follows that Tk is a monotonous 
function of D/fi. By choosing ttuq = ^InD in (|T1T|) . one arrives at the equations yielding the free 
energy of a Kondo impurity, 

Inb(x) = -e^ + /3/i/2+ [A:*ln5S](x) - [A;*ln^](x + i7r-ie) (26) 
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and *B = Thus, the functions 55, !B are decoupled from <t. This decoupUng rests on H24|) . 

which is vahd for [Inc'] 3> 1 and Ap <C (lnL')/7r = uq at fixed n/D. The constant Tk which 

appears in (|27j) is the scahng invariant expected from Poor Man's scahng in the Kondo regime. This 
becomes evident by writing 

Tk = Tk/D = TKexp[-Tr\uo\], (28) 

where from (|25|) . k oc D/T. The freedom in the choice of the sign of uq corresponds to the appearance 
of antiferromagnetic exchange with the right or left movers in ©, following the sign of uq. By 
comparing (|25j) with (|19j) . one sees that 7 = DTk- 

If ^ 1 but still finite, from the NLIE an approximate expression for Ind is obtained, 

,/3/i/2-/3<}!>cM/2+/3m , „-/3V2-/3</>c(«)/2+/3m^ / ( m(oo)«^{°o) 



where is defined in (gSJ in appendixEland !B(°°) = l + exp(/3/i + /3/i/2), = 1 + exp (/?//- /3/i/2)- 

Then from ()48|1 in appendix 

T /■°° ine:(u') , , , 

uncharge = / 7 TTTT^"'' (3^) 

which is the free energy of an effectively free impurity with a certain density of states, each state 
permitting no and single occupation. The approach to the Kondo susceptibility is depicted in fig. [21 
for constant = 2/5. 

Apart from the susceptibility, it is instructive to consider the specific heat at constant ^, shown 
in fig. 131 The approach of the Kondo limit reveals two different scaling regimes. In a first step, for 
Uq < Ac, spin and charge excitations are not decoupled, but occur on the same energy scale. This 
is shown by the single peak in fig. |31for low enough D. With increasing D, the shift uq increases; 
once Uq > Af, the charge excitations are scaled to higher and higher energies ~ D, whereas the new 
universal scale Tk (only dependent on the ratio D/fj,) appears as characteristic scale for the spin 
excitations. For uq S> A^, spin- and charge excitations are energetically well separated as described 
by (|23), JSni). 

3.4 Phase diagram 

The results of the previous sections 13.1113.31 are summarized in figure |H1 For low temperatures, the 
three different impurity regimes are sketched in the fi — D-plane with uq = InD/vr. The LM-regime is 
distinguished from the other two by the straight line fj, = Hc, for fj^c, see (fT^ . The SC- and MV-regimes 
are separated by the condition \uo\ = InA^. From (|T7|l . 

fi = 4:D(a + l)/a(l ^ In^ D ) (31) 

for |Ac| <^ 1 and D > 1. Note that in (|31|) . fi ^ fic for D — > 1. In the other extreme, at |Ac| ^ 1, the 
driving terms in the linearized equations are expanded appropriately, yielding 



A^ = D{a + l)a/n 
H = D-K'^{a + l)a/\T? D, ioT uo = K. (32) 
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2 3 4 5 6 7 

InT 



Figure 1: Effective Curie constant of the impurity at ^ = 300 for different D. The curves with 
D = 1,30 exhibit LM-behavior, the one with D = 300 corresponds to the MV-regime. In the other 
cases shown, signatures of the SC-regime are discernible. At high temperatures, the asymptotic values 
1/6 and 1/4 are j 




Figure 2: Impurity susceptibility Xi and T ■ Xi at constant D/^ = 2/5 illustrating the crossover from 
Anderson to Kondo behavior. At high temperatures, the asymptotic values 1/4 and 1/6 are indicated. 
For the ratio D/^ = 2/5, the zeroes of Inc are found numerically to be zizAp = 0.893. So uq < Ap 
corresponds to D < 16. Note the different scaling behaviors for uq ^ Ap. 
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Figure 3: Impurity specific heat at constant D/^ = 2/5 in the crossover region from Anderson to 
Kondo behavior. The inset shows the linear T coefficient at low temperatures, the Kondo value 
being 7r/3. The peak corresponding to charge excitations shifts to higher temperatures. For the ratio 
D/ ^ = 2/5, the zeroes of Inc are found numerically to be ztAi? = 0.893. So uq < Ap corresponds to 
D < 16. Note the different scaling behaviors for uq ^ Ap. 




Figure 4: Impurity specific heat at constant D/iJ, = 2/5 in the crossover region from Anderson to 
Kondo behavior, with the temperature scale T/D. Also shown is the specific heat derived from (|29|) . 
(|5n)). that is an effectively free three-state impurity, with parameters D = 3200, n = 8000. 
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Figure 5: Phase diagram showing the three different regimes of the impurity (diagonally hatched: LM, 
only vertically hatched: SC, vertically and horizontally hatched: MV) at low temperatures. The axes 
begin at D = 1, //q := ^J'c{D = 1) = 4(a + l)/a. The line between MV and SC is the center of a 
transition region between these two regimes, determined by the condition uq = InD/n. The impurity 
behaviour along the horizontal (diagonal dashed) line at fi =const. {f^/D =const.) is illustrated in fig. 
in pilH) . In the host, only two phases are present, LM (diagonally hatched) and MV (vertically hatched, 
that is MV and SC of the impurity). In the latter region, the host behaves like a Luttinger-liquid, 
whereas it is completely filled for fi > fic- 



Eq. shows that the SC-regime can always be obtained for n < fic, fJ'/D fixed in the limit I? — > oo. 
Also shown in the phase diagram El are the regions illustrated by figures HEl 

The SC regime can only be realized for the impurity. As far as the host is concerned, the choices 
^ > muc (p < fJ-c) lead to a completely (not completely) filled band. In the second case, the host 
behaves as a Luttinger liquid at low temperatures, cf. appendix IbI 

3.5 Luttinger-liquid-behaviour of the host 

Since uq = in the host, the only temperature scale is the bandwidth D, and all thermodynamic 
functions depend only on T/D at constant D/ii. In the LM-regime of the impurity, the band is 
completely filled, allowing for no excitations of the host. Away from complete filling, the host shows 
Luttinger-liquid-behavior for f3D S> 1, as shown in appendix IbI 

A = _r^=!(i + 1) 

6 \Vc VsJ 

with distinct spin- and charge velocities Vg and Vc- This corresponds to a linearized energy-momentum 
dispersion in 0, the host interactions still being finite. In the Kondo limit, i.e. ^ oo at constant 
D/iJ,, the Luttinger-liquid regime of the host is valid for all temperatures, including, as a special case, 
the free fermion limit at q — > oo. Thus the model presented here exhibits Kondo-like behavior even 
in the presence of an interacting bulk. In [12, 13] arguments of CFT have been presented in order 
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to decide whether generic singlechannel Kondo-like behavior is possible in an interacting host or not. 
The answer to this was found to depend on the host's properties. Further studies of our model with 
techniques of conformal field theory are expected to yield additional criteria to answer that question. 

4 Conclusion 

We have studied an asymmetric Anderson model in which a three-state impurity permits charge and 
spin exchange with an interacting bulk on a lattice. For this model, we distinguished a mixed valence, 
a strong coupling and a local moment regime with the help of scaling arguments. This distinction was 
made quantitative by the exact calculation of the free energy contribution of both the impurity and 
of the host. On the basis of the exact solution, the three different regimes have been studied both 
analytically and numerically. Analytically, a new universal energy scale for excitations was deduced. 
Its gradual emergence was illustrated by numerical calulations. The numerical solution illustrated the 
emergence of a new energy scale, namely the Kondo temperature Tk, in the strong coupling limit. 
In the Kondo limit, this scale is separated by an increasingly large energy from the temperature on 
which charge fluctuations are excited. Our model can be generalized to a four-state impurity, which 
also permits double occupation. The derivation of the Hamiltonian and its analysis, especially the 
comparison with the finite- [/ Anderson model, are projects for future work. 



A Appendix A 



In this appendix we sketch the construction of the gl(2|l)-symmetric transfer matrix, from which the 
Hamiltonian is obtained, and the derivation of the NLIE (fTn|) - (|12|) . 

Let be the module of the 4-dimensional irrep of gl(2|l), labeled by the parameter a. For 
a = 0, this is the three-dimensional fundamental irrep. Following [16,25], we define matrices 

gEnd(^l/('^) «) where the first (second) space is referred to as "auxiliary" ("quantum") 

space. These matrices satisfy the graded Yang-Baxter-Equation 



/?,7 



13', -y' 



5,7' 



5', 7" 



5', 13' 



i?("'"')(^z) i?(""'"')(^;) " R^-"'''\v-u) (-l)m+P[S'])Pm 

13", y" 



5,(3 



5', 13' 



r('^"^»)(v-u) ' i?(""'"')(7;) ' i?("'"')(n) 



5', 7 



5", 7' 



(33) 



Here a grading is assigned to the basis vectors through the parity function p, namely p[l] = p[4] 
0; p[2] = p[3] = 1. The matrix i?^"'" ^ is written as a sum of projectors. 



1 + ^P,("'"') 



_J_p{a,a') 

p{v) 3 



where j5(u) = l/2 + t>/(a + Q' + 2), q{v) = — l/2 + ?;/(a + a') and Pi, P3 are projectors from ) 
onto gl(2|l) modules with highest weights (0, 0|a + a') and (— 1,— l|a + a' + 2) respectively. They 
are given explicitly in [16,25]. Note that i?^"'" \v) R^°' '")(—?;) = 1. Furthermore, consider the set of 
matrices 



(u) 



7,5 



(-1) 



p[5]{p[^]+p[a]) 



-u 



1,0 



(34) 
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which also satisfy a YBE. From R and R, monodromy matrices T and T are constructed by multiplying 
(L + l)-many i?-matrices in auxiliary space, 

T{u) = d:x^^-\u)d:x-i^{u)...R^:-^^^^^ (35) 

T(n) = Rir'\-u)R^:X'^^\-u)...R^^^^^^ 

where indices a {u = 0, . . . , L) refer to the auxihary (quantum) spaces. Note that the quantum space 
of the last factor (impurity) carries a representation different from that in the host, ai ^ ah. For the 
model (pQ), one chooses = 0, ah = a. A shift by iuq of the spectral parameter is allowed on the 
impurity site. By taking the supertrace over the auxiliary spaces of T, T, one obtains the transfer 
matrices r, r, giving rise to the Hamiltonian H, 

t{u) = straT(n) , r(?i) = straT(n) (36) 
In [rr] (n) = In [rr] (0) + u [r-i(0)r'(u) + r~Ho)r'(u)]^^Q +0{u^) . (37) 

^ V ^ 

= : const. if 

One is free to choose the multiplicative constant in the definition of H by scaling u. Here, we take 
-D(a/i + l) as a common prefactor, where D is a bandwidth parameter. External fields ^, h are included 
by twisted boundary conditions of the quantum spaces, for details on the calculation of see [9]. 

Consider the quantum transfer matrix rfj^ pertaining to the impurity/host sites, [9], with Trotter 
number A^. The free energy contributions of the impurity and host are given by the largest eigenvalue 
A™^ of the QTM in the limit of infinite Trotter number, 

/, = - hm 1 In Ar^(no) , fh = - lim \ In A)r^(O) . (38) 

All eigenvalues of t^^^ can be found by ABA-techniques similar to [8, 15, 16], so that we do not repeat 
the calculation here. The result is 

A.'«(„) = ^.(„)*(!i±i±Me'* 

qi [v + 2) 



qi {v + iah + I) g2(^ + i) ^f3h/2 ^ gi {v + + ^) g2(-^ - i) 



qi {v + iQfe + 2) 
qi {v - 2) 

The vacuum expectation values are 



[V — 1- 

{v - - i^) cb+{v + i + i^ + if) (/>_ {v + i^ + if) 
Av-i^ + if) ^. {v-i^+if) 



(l)+{v + i + i^ + if) (/._ {v + i^ + if) 

Mv) = M-v)L,^-a,,^-l ' 
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where 4'±{'^) ■= (v i iu)^^'^ and u := —D(ah + 1)P/N. The g-functions are polynomials with zeroes 
at the BA-numbers, Quiv) '■= Y\k=i ~ '^i'''*)' where u = 1,2. The two sets of BA-numbers |^fc^^|) 
I'^fc^^} determined by the analyticity of A'p\v), 



q2 [v'^^ - i) qi (v^k^ 



The eigenvalues of A^^^^^ of the host matrix are given by an expression similar to ()39p. with Qj = a/j. 

The free energy follows from the largest eigenvalue, determined by Mi, = N/2 many BA-roots 
in both sets i/ = 1,2. For this case, lnA™'^^(f) has been calculated in [20,27] from appropriately 
chosen auxiliary functions, obeying a closed set of non-linear integral equations (NLIE). Within this 
approach, the Trotter limit — > oo can be carried out analytically. It turns out that lnAf^^^(v) can 
be calculated analogously to lnA^^^{v) with the result: 

lim lnAf'^^(z;) = r]{v) + [(* In ^]{v) + [C * ln'^]{v) + + In (t\{v) , (40) 

A''— ►oo 

where the auxiliary functions 03 = l + b,53 = l + b and <L = 1 + c obey the NLIE 

lnb(t;) = Mv + iS)-[h*ln'^]{v + 2i6)-[ki,*ln(!:]{v + i6)+p{fi + h/2) (41) 
lnb(?;) = (p^{v-i6)-[k-^*ln^]{v-2i6)-[k-^*ln^]{v-i6)+/3{^x-h/2) (42) 
lnc(?;) = ^c{v)-[kb*ln'^]{v + id)-[k-^*ln^]{v-id)-[k,*ln(i\{v) + 2p^i. (43) 

Here, the driving terms 

{v ^ \ahl2){v - \ahl2 - I) 
(/)c = (/'b + '^fa (45) 

and convolutions [/ * := /(w — w)g{w)dw with integration kernels 



1 

2ttv{v — i) 

{ah + l)iai/2 + ah/2 + 1) 



kb(v) — 7^ZZ77. IT ' ~ ' — ^b + ^t, 



a,^a, v^ + {a,/2 + ah/2 + ir 



Dpiah + 1) 

have been defined. The external fields enter through the asymptotic values of the auxiliary functions, 

g/3(M+/i/2) _ e/3(M-h/2) 

b(±oo) = 37 — 7-17^ ,b(±oo) = „, , , (46) 
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Prom H38|) and (|4U() . the free energy contribution of the impurity results into 

_ J -r7?M-r[C*hi5S](no)-T[C*hi^](no)-T[(C + C)*ln(2:](no), / 

" I T In c (no) + 4D(a;, + 1)^^-2^, «i = 0. ^^^^ 

The free energy of the host per lattice site is also given by eq. (|48|). with ai ^ at, 

fh = fi\ai^ah;uo=0 ■ (^9) 

Note that, apart from the //-term in ()4U() . the parameter a/^ enters explicitly only in the NLIE 
(|1T|) - (|43|) . whereas is found only in the free energy equation (|i8|) . This had to be expected: The 
quantum spaces of t-^^ are modules of the a/i-representation of gl(2|l). The auxiliary space, over 
which the trace has to be taken to get t^^\ carries the Oj-representation. Especially, the auxiliary 
functions depend only on ah, and thus are identical for the host and impurity eigenvalues. However, 
the functional dependence of the eigenvalues on the auxiliary functions is different for the host and 
the impurity. 

The NLIE (UHIl-lIll) are obtained from by setting 6 = 1/2, ah = a, at = and by making 

the substitution 



[ki,,c * In e:] (v) = [ki,,c * In c] (v) + 



c 



Manipulations in Fourier space yield the result (fT?H) - (|T2|) . For Oj ^ 0, eq. (|i8|) yields the free energy 
of a four-site impurity, where the charge fluctuation terms include transitions from and to the doubly 
occupied state. 



B Appendix B 

We calculate the contribution 0{T^) to fh, eq. To this end the NLIE are linearized for 

(3D,(3fj,,P » 1, following a similar procedure applied to the homogeneous tJ-model in [19]. Since 
no = 0, the form H4H) - (|43|) is best suited for our analysis. First observe that for h sufficiently large, 
b = O (e~^'^) and therefore *B = 1 up to exponentially small corrections. Secondly, both Inb, Inc have 
two zeroes, lnb(±A[,) = and lnc(±Ac) = 0. Due to the driving terms in (|IT|) - (|1S)) of order 0{P), 
the slopes of the auxiliary functions in their zeroes behave as | In'j, |.(ibAt,,c)| ~ P, so that there is a 
sharp crossover for /? — > co in the region around ±A(,,c- Thus we define the scaling functions 8^^^ by 
j3eb := Inb, f3£c '■= Inc. Then 

lnQ3(n) = Pe,{v)e {v^ - A^) + ^^^^J^'^^^^^ {6{v - A,) + 6{v + A^)) , (50) 

and similarly for \n<L, [19,27]. Here, 9{v) denotes the Heaviside function. Higher order terms can be 
neglected if the slope is sufficiently steep. This approximation is inserted into H4H) - ()43|) with 6 = 1/2, 
such that the functions are real-valued on the real axis. With Kb,c{'^) •= ^fa,c(f^ + i/^)? this gives 



^biv) = £o,biv) - / Ki,{v - w)£c{w)dw (51a) 

J\w\>Ac 



Eciv) = eo^c{v) — / Hib{v — w)£ii{w) dw — / K^{v — w)ei:{w) dw , (51b) 
'\w\>At, J\w\>Ac 
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where we have defined the dressed energy terms 

, , (hbiv + i/2) h vr^ r , k \ /am 

eo,{v) := + ^ + 2 - l2(i^K{K)\ ^^^^^ " + ^^^^ + ^^^^ 



7r2 



K,(u-A,) + K,(v+Ac)] . 



12/?2|£UAcyi 
Eqs. (|51a|) . ()51b|) are written compactly as 

1 = ( ) ' (52) 

Kb 1 + y V y V '^o.c y 

where the convolutions are to be taken in the appropriate hmits. Define two further functions ^b,c by 

1 «b = ( h 

/^b 1 + y ) \ fc 

with the same integration hmits as in eq. 1)52(1 and 

/b : 



It follows that 





Oh 




«^b, 


Oh 


= 




-ah ' 


a/i 7^ 






ah = 0. 



/c := <^ / (53) 



/ eb(i')/b(^^) + / ec{v)f^{v)dv 



'j\>Ab Jlvl 

Cb{v)eo^b{v){v)dv+ / Uv)eo,c{v) dv . (54) 

?)|>Ae J\v\>Ac 

Writing the impurity contribution to the free energy for uq = 0, the second case in ()48)). in terms of 
the scaling functions (|5(J() . and using l|54p . one obtains 



V ^y i|i.|>Ai, J\i 

-T [ 4>^{v + \/2)i^{v)dv-T f 4>,{v)i 

J\v\>ki-, Jv>Ac 



\v\>Ab -'|t'|>Ac 



'v) dv 



This result reflects that at low temperatures, the contributions of charge and spin excitations can be 
separated, but occur on the same temperature scale. The charge and spin velocities are given by 

1 _ gc(Ac) 1 _ gb(Ab) 

For the explicit calculation of Vc,s in certain limiting cases in the range of A[,,c) one may use Wiener- 
Hopf techniques; this has been done for Oi = Oh = [28] and a/i = 0, Qj / [8]. The calculations of 
this appendix also apply to the D(3 ^ 1-behavior of fi, if uq = 0. This demonstrates that the shift 
uo ^ liiD induces a new energy scale which energetically separates charge from spin excitations. 
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